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We present a multipartite entanglement measure for iV-qudit pure states, using the 
norm of the correlation tensor which occurs in the Bloch representation of the state. We 
compute this measure for important class of iV-qutrit pure states, namely general GHZ 
states. We prove that this measure possesses almost all the properties expected of a 
good entanglement measure, including monotonicity. Finally, we extend this measure to 
./V-qudit mixed states via convex roof construction and establish its various properties, 
including its monotonicity. 



1 Introduction 



Entanglement has proved to be a vital physical resource for various kinds of quantum 
information processing, including quantum state teleportation [1,2], cryptographic key 
distribution [3], classical communication over quantum channels [4,5,6], quantum error 
correction [7], quantum computational speedups [8] and distributed computation [9,10]. 
Further, entanglement is expected to play a crucial role in the many particle phenomena 
such as quantum phase transitions, transfer of information across a spin chain [11,12] etc. 
Therefore, quantification of entanglement of multipartite quantum states is fundamental 
to the whole field of quantum information and in general, to the physics of multicomponent 
quantum systems. Whereas the entanglement in pure bipartite states is well understood, 
classification of multipartite pure states and mixed states, according to the degree and 
character of their entanglement is still a matter of intense research [13,14.15]. Principal 
achievements are in the setting of bipartite systems. Among these, one highlights Woot- 
ter's formula for the entanglement of formation of two qubit mixed states [16], which still 
awaits a viable generalization to multiqubit case. Others include corresponding results 
for highly symmetric states [17,18,19]. 

Interest in multi-dimensional entangled states comes from the foundations of quantum 
mechanics as well as the development of new protocols in quantum communication. For 
example, it has been shown that maximally entangled states of two quantum systems in 
a high dimensional Hilbert space, qudits, violate local realism stronger than qubits, and 
entangled qudits are less affected by noise than entangled qubits [20,21]. In quantum 
cryptography [22], the use of entangled qutrits [22,23] or qudits [24,25] instead of qubits 
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is more secure against eavesdropping attacks. Furthermore, the protocols for quantum 
teleportation or for quantum cryptography work best with maximally entangled states. 
These facts motivate the development of techniques to generate entangled states among 
quantum systems in a higher dimensional Hilbert space with good entanglement charac- 
teristics. Technical developments in this direction have been made. For example, four 
polarized entangled photons have been used to form two entangled qutrits [26]. Entangled 
qutrits with two photons using an unbalanced 3-arm fiber optic interferometer or photonic 
orbital angular momentum have been demonstrated [27,28]. Time-bin entangled qudits 
of up to 11 dimensions from pump pulses generated by a mode- locked laser have also 
been reported [29]. In short, quantifying the entanglement measure of a qudit system is 
of physical interest. The issue of entanglement in multipartite and higher dimensional 
states is far more complex. Notable achievements in this area include applications of the 
relative entropy [30], negativity [31] Schimidt measure [32] and the global entanglement 
measure proposed by Meyer and Wallach [33]. 

A measure of entanglement is a function on the space of states of a multipartite system, 
which is invariant on individual parts. Thus a complete characterization of entanglement is 
the characterization of all such functions. Under the most general local operations assisted 
by classical communication (LOCC), entanglement is expected to decrease. A measure 
of entanglement which decreases under LOCC is called an entanglement monotone. On 
bipartite pure states the sums of the k smallest eigenvalues of the reduced density matrix 
are entanglement monotones. However, the number of independent invariants (i.e. the 
entanglement measures) increases exponentially as the number of particles N increases 
and complete characterization rapidly becomes impractical. A pragmatic approach would 
be to seek a measure which is defined for any number of particles (scalable), which is 
easily calculated and which provides physically relevant information or equivalently, which 
passes the tests expected of a good entanglement measure [13,14]. 

In this paper, we present a global entanglement measure for iV-qudit pure states which 
is scalable, which passes most of the tests expected of a good measure and whose value 
for a given system can be determined experimentally, without having a detailed prior 
knowledge of the state of the system. The measure is based on the Bloch representation 
of multipartite quantum states [34] . 

The paper is organized as follows. In section 2 we give the Bloch representation of a 
./V-qudit quantum state and define our measure E?. In section 3 we compute E*r for an 
important class of iV-qudit states, namely, GHZ states. In section 4 we prove various 
properties of Ej-, including its monotonicity, expected of a good entanglement measure. 
In section 5 we extend E? to iV-qudit mixed states via convex roof and establish its 
monotonicity. Finally, we summarize in section 6. 
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2 Bloch representation of a iV-partite quantum state 



Bloch representation [35,36,37,38,39] of a density operator acting on the Hilbert space of 
a d-level quantum system C d is given by [40] 

i 

Eq.(l) is the expansion of p in the Hilbert-Schmidt basis = 1, 2, . . . , d 2 — 1} 

where Aj are the traceless hermitian generators of SU(d) satisfying Tr(XiXj) = 25^ and 
are characterized by the structure constants of the corresponding Lie algebra, fijk,gijk 
which are, respectively, completely antisymmetric and completely symmetric. 

2 

— ~^ i jl< 1 + ifijk^k + Qijk^k (2) 

s = (si, S2, • • • , Scp-i) in Eq.(l) are the vectors in IR 6 ' 2 " 1 , constrained by the positive 
semidefiniteness of p, called Bloch vectors [38]. The set of all Bloch vectors that consti- 
tute a density operator is known as the Bloch vector space _B(IR d2_1 ). The problem of 
determining B(R. d ~ 1 ) where d > 3 is still open [36,37]. However, for pure states (p = p 2 ) 
the following relations hold. 



IH| 2 = \j d<yd 2 SiSjgtjk = (d - 2)s k (3) 

where ||.||2 is the Euclidean norm in M d2_1 . 

It is Known [41,42] that B(M d2 ~ 1 ) is a subset of the ball D R (M. d2 ' 1 ) of radius R = 

yj ^ d( - d ~^ ) which is the minimum ball containing it, and that the ball D r (]R d2_1 ) of radius 
r = is included in B(R d2 ~ 1 ). That is, 

D r (R d ' 2 ^ 1 ) C B(R d2 ^) C D^R*- 1 ) (4) 

In order to give the Bloch representation of a density operator acting on the Hilbert 
space C d <E> C d <g> • • • ® C d of a iV-qudit quantum system, we introduce following notation. 
We use k, ki (i — 1, 2, • • • ) to denote a qudit chosen from iV qudits, so that k, ki (i = 
1, 2, • • • ) take values in the set M = {1, 2, • • • , N}. The variables a k or a kl for a given 
k or hi span the set of generators of SU(d) group (Eqs.(l) and (2)) for the kth or fcjth 
qudit, namely the set {Ai, A 2 , • • • , \d 2 -i} f° r the k^th qudit. For two qudits ki and k 2 we 
define 

^ = (h ® Id ® • • • <8> A afci <g> J d <g> • • • (8) J d ) 
A£ 2) = (/ d ® / d ® • • • ® A Qfc2 <g> J d (8) • • • (8) J d ) 
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x ££ x %£ = (ld®ld®---® K kl ®h®---® K k2 ®h® h) (5) 

where A afci and A Qfe2 occur at the A^th and A^th places (corresponding to fcith and 
/c 2 th qudits respectively) in the tensor product and are the a^th and «fc 2 th generators of 
SU(d), a kl = 1, 2, . . . , d 2 — 1 and a k2 = 1, 2, . . . , d 2 — 1 respectively. Then we can write 

' = ^W + EE«-*S? + E E <«w£Xt l + ---+ 

k&M a k {ki,k 2 } afe 1 Ofe 2 

V V t A (fel) A (fe2) ---A (feM) +---+ V t A (1) A (2) ---A (JV) I 

{fci,fe 2 ,-- ,&m} a fc 1 a fe 2 '" a fcM a 1 a 2 -a N 

(6) 

where is a Bloch vector corresponding to Arth qudit, = [s afc ]a fc =i "which is a 
tensor of order one defined by 

*«* = %rr\p\<£] = ^Tr[ Pk X ak ], (7a) 

where p k is the reduced density matrix for the fcth qudit. Here {ki,k 2 ,--- , A;^}, 2 < 
M < N, is a subset of A/" and can be chosen in (^) ways, contributing (j^j terms in the 
sum E^/ fcl fca ... fcM } in Eq.(6), each containing a tensor of order M. The total number of 
terms in the Bloch representation of p is 2 N . We denote the tensors occurring in the sum 

E {klM ,..., kMV (2 < M < N) by T<*^-' fc "> = [t QfciQfe2 ... afeM ] which are defined by 

+ = — Tr\o\^ kl) \ {k2) ■ ■ ■ A (fcM) l 

L a kl a k2 ...a kM 2 M lF ak l ak 2 a *M J 

= ^M Tr Kfc2...feM( A « fcl ® V 2 ® • • • ® V M )] (7b) 

where pk x k^...k M is the reduced density matrix for the subsystem {k\k 2 . . . ku\- Each 
of the (^) tensors of order M, occurring in the Bloch representation of p, contains all 
information about entanglement of the corresponding set of M subsystems. All informa- 
tion on the entanglement contained in p is coded in the tensors occurring in the Bloch 
representation of p. The tensor in last term in Eq. (6), we call it T^ N \ contains all the 
information of genuine iV-partite entanglement. 

The following we give the generators of SU(3) in the |1), |2), |3) basis for later use. 
A 1 = |l)(2| + |2)(l| 

A 2 = -i(|l)(2|-|2><l|) 
A 3 = |l)(l|-|2)(2| 
A 4 =|l)(3| + |3)(l| 
A 5 = -.(|l)(3|-|3)(l|) 
A 6 = |2)(3| + |3)(2| 
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A 7 = -i(|2)<3|-|3)<2|) 

A 8 = ^(|l)(l| + |2)(2|-2|3)(3|) 

The action of these generators on the basis states {|1), |2), |3)} is given by 

Aijl> = |2), Ai|2> = |1), Ai|3> = 
A 2 |l) = i|2), A 2 |2) = A 2 |3) = 

A 3 |l) = |l), A 3 |2) = -|2>, A 3 |3) = 
A 4 |l) = |3), A 4 |2> = 0, A 4 |3) = |l) 
A 5 |l)=i|3), Ax|2> = 0, Ai|3> = — i|l> 
A 6 |1) = 0, A 6 |2) = |3), A 6 |3) = |2) 
A 7 |1) = 0, A 7 |2)=i|3), A 7 |3) = -i|2) 
A 8 |l) = ^|l), A 8 |2) = ^|2), A 8 |3) = -^|3) 

We will use these equations below (see the next section and proof of proposition 7). 
We propose the following measure for a iV-qudit pure state entanglement 

^>)HI^I-(^^) W2) (8) 

where is given by Eq.(7b) for (M = N) in Bloch representation of p = \ip)(tp\. 
The norm of the tensor appearing in definition (8) is the Hilbert-Schmidt (Euclidean) 
norm ||TW|| 2 = (T^,^) = Z) aiaa ... Qjv fa iaa ... ajv - Throughout this paper, by norm, 
we mean the Hilbert-Schmidt (Euclidean) norm. 



3 Testing out Et(\iP)) for qutrit pure states 

Before proving various properties of Er(\ip)) we calculate it for 3- and 2-qutrit pure states 
employed to study multipartite entanglement in the literature. The most important 3- 
qutrit state is 

|V>> = a|lll) + (3\222) + 7 |333), a 2 + f3 2 + 7 2 = 1 (9) 

where |2), |3)} is the computational basis for each qutrit. This state and its 2-qutrit 
restriction 

|V>) = a|ll> + (3\22) + 7 |33), a 2 + (3 2 + 7 2 = 1 (10) 

are used in various ways to understand qutrit entanglement. It is shown to break Clauser- 
Horn-Bell type of inequality for three qutrits [27]. It is used in protocols for quantum key 
distribution based on encoding in qutrit systems [23]. The 2-qutrit version is experimen- 
tally prepared [28,29]. Below we get a closed form expression for Ej- for 3-qutrit state, 
Eq.(9) and we relate E? for the 2-qutrit state, Eq.(10), with 2-qutrit concurrence [43]. 
For the state in equation (9), the general element of T N is given by 

3 N 

thi 2 -i N = t^&IKK ■ ■ -KM ik = • ,8, k — 1,2,3. 
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From the action of A operators on the basis states as shown in the previous section, there 
are only 20 nonzero elements of T N . We have t 888 = -^[a 2 + (3 2 + (— 2) 3 7 2 ]. t m = 2a/?, 
^333 = a 2 — j3 2 , £444 = 2cry, £666 = 2/3 7 . Other non-zero elements correspond to two A2S 
and Ai; two A5S and A4 and two A7S and Xq. Further, t 833 = t 383 = t 338 = [a 2 + [3 2 ] and 

^838 — ^388 = ^883 = 3 I® 2 — P 2 ]- Note that for a = (3 = 7 = we have only 16 non-zero 
elements. Thus the expression for Eq- becomes 



Mm) = \\T {N) \\ - 3 3/2 = ^[27(a 2 /5 2 + «V + PV) + -M 2 - ff + ^(a 2 + p 2 ) 2 

A 4 lb 

+ ^(a 2 + /3 2 - 8 7 2 ) 2 ] 1/2 -v / 27. (11) 
The value of E T for a = p = 7 = ^ is 3.0196. 



Y:33 
Z: 3.016 




The N-qutrit state corresponding to Eq.(9) for iV-qutrits is 

|V>) =a|lll---l) +/?|222---2) + 7 |333---3) a 2 + (3 2 + 7 2 = 1 (12) 
Then any element of namely 



3iV 



U li2 ...i N = -^wi^lKK ■ ■ • A ijv |-0>, ifc = 1,2, • •- ,8, A; = 1, • • • ,N (13) 
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which has 2k A 2 s and (N — 2k) Ais k — 0, 1, • • • , [yj, where |_^J is the greatest integer 
less than or equal to x, is equal to (— l) k 2a[3. If the element has 2k A 5 s and (N — 2k) A 4 s 
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k — 0, 1, • • • , LyJ then the value is (— l) fc 2cry. Similarly, if the element has 2k A 7 s and 
(N — 2k) A 6 s then the value is (-l)^/^. Further, if the element has 2/cA 3 and (N — 2k)\ 8 
with fc = 1, 2, • • • , [f J then its value is (-^) N - 2k (a 2 + (3 2 ). If the element has (2k + 1)A 3 
and (N — 2k — 1)A 8 with fc = 0, 1, 2, • • • , [f J - 1 then the value is (^) Ar_2fe_1 (Q; 2 - (3 2 ). 
Now t 88 ... 8 = (-^) N (a 2 + (3 2 + (-2) Ar 7 2 ). All other elements are zero. Therefore, 

LfJ LfJ 

\\ TN h 2 = ( ^) 2 [ E ( 2 1) ( 4 ^ 2 + 4 « v + 4/5 V) + E (*) dr 2 v + /? 2 ) 2 

fc=0 ^ ' fc=l ^ ' 

LfJ- 1 

+ E L+O^ -2 ^^ 

fc=0 ^ ' 

£^(|^)) = ||T JV ||-3 JV/2 . (14) 

We can check that Eq- for in Eq.(14) for three qutrits is the same as Eq- for 
in Eq.(ll) for N qutrits with N = 3. When N = 2, i.e. two qutrits with given by 
Eq.(10) we get 

Erm = (l)[S(a^ + aW + /?V) + (1 - 7 2 ) 2 + (jj)(« 2 - + (±)(1 + 3 7 2 ) 2 ] 1/2 - 3. 
The concurrence for the state for two qutrits in Eq.(10)is [43] 

C{\4))) = v/W + aV + ^V). 
Therefore, for this two qutrit state, 

E T = ( 9 -)[2C\m + (1 - 7 2 ) 2 + (jj)(a 2 - /? 2 ) 2 + (±)(1 + 3 7 2 ) 2 f/ 2 - 3. (15) 
We can check that whenever the concurrence C = 0, Eq- = 0. 

4 Properties of Et{\^)) 

To be a valid entanglement measure, Eq-(\ip)) must have the following properties [44,45]. 

(a) (i) Positivity : Eq-(\ip)) > for all iV-qudit pure state \ip). (ii) Discriminance: 
Et(\iP)) = if and only if \tp) is separable (product) state. 

(b) LU invariance : Eq-(\ip)) is invariant under local unitary operations. 

(c) Monotonicity : local operators and classical communication (LOCC) do not in- 
crease the expectation value of Eq-(\ip))- 

We prove the above properties for Eri\ip)). We also prove the following additional 
properties for Eq-(\?p)). 
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(d) continuity ||(|^)^| - |0)(0|)|| ^ ^ E{\$)) - E(\<j>)) - 0. 

(e) superadditivity E T (\ijj) <g> |</>)) > £r(|^)) + Et(\4>))- 

We need the following result which we have proved in [34] . 

Proposition : A pure iV-partite quantum state is fully separable (product state) if 
and only if 

rw = s« o s ( 2 ) o . . . o s w = o£U* (fc) , (16) 

where s( fc ) is the Bloch vector of kth subsystem reduced density matrix. 

The symbol o stands for the outer product of vectors defined as follows. 

Let uW, u( 2 \ . . . , u( M ) be vectors in R d ? _1 , R^' 1 , ■ ■ ■ , R^m" 1 . The outer product o 
u^ 2 ^ o • • • o u^ M ^ is a tensor of order M, (M-way array), defined by 

f . . . _ u (i) u (2) m. i < < fp _ i._io... m 

Proposition 1 : Let be a A^-qudit pure state. Then, HT^^H = (d(d — 1)/2)( N / 2 ^ 
if and only if is a separable (product) state. 

Proof. By proposition 0, \ip) is separable (product) if and only if 

As shown in [46,47 ], 

{OLyS (k \OLrS (k) ) = ^M k \s { % (17) 
where (, ) denotes the scaler product. This immediately gives, for qudits, 

l| T W||2 = (T w tW) = nf =1 ( s « s ( fc )) = n fe || s «|| 2 = ( ^_il2y 



Proposition 1 immediately gives 

Proposition 2 Let \ip) be a A^-qudit pure state. Then Er(\ip)) = if and only if 
is a product state. 

Proposition 3 : Let be a N-qudit pure state. Then \\T^ N} \\ > (d(d - 1)/2)W 2 ). 

If |^) is not a product of A" single qudit states (i.e. is not A^-separable) then it is 
N — k separable k — 2, 3, • • • ,N — 1. Viewing the A^-qudit system as a system comprising 
N — k qudits, each with Hilbert space of of dimension d and k entangled qudits with the 
Hilbert space of dimension d k , we can apply proposition to this separable system of 
N-k + 1 parts in the state We get = (s (1) ) o (s^) o • • • (s {N ' k) ) o ( s ( N - k+r >) 

This implies, as in proposition 1, via Eq. (17) and Eq. (3) that 

ll^ll 2 = nf_f +1 ||( S «)|| 2 = d ^- l ) (M_Af- k > ( 4 = d*). (is) 

If k = N we attach an ancilla qudit in an arbitrary state \<f>) and apply proposition to 
A^ + 1 qudit system in the state <g> \(f>) where \ip) is the A^-qudit entangled state. This 
result, combined with proposition 1, completes the proof. 
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Proposition 3 immediately gives 



Proposition 4 ■ E T (\ijj)) > 0. 

We now prove that E*r(\ijj)) is nonincreasing under local operations and classical com- 
munication. Any such local action can be decomposed into four basic kinds of operations 
[48] (i) appending an ancillary system not entangled with the state of original system, (ii) 
performing a local unitary transformation, (iii) performing measurements and (iv) throw- 
ing away, i.e. tracing out, part of the system. It is clear that appending ancilla cannot 
change \\T^\\. We prove that Er(\ip)) does not increase under the remaining three local 
operations. 

Proposition 5 : Let Ui (i — 1, 2, • • • , N) be a local unitary operator acting on the 
Hilbert space of ith qudit . 
Let 

p = (CiWCCiC?) (is) 

for density operators p and p' acting on Tt = ®f =1 H^ and let T^ N "> and T'^ denote the 
N partite correlation tensors for p and p' respectively. Then, 
\\T' {N) \ \ = \\TW\\, so that E T (p) = E T (p') 

Proof. Let U denote a one qudit unitary operator, then, there exists an orthogonal 
matrix [0 Q/ g] acting on R d such that UXaW = YlpOaP^P 

where [O a p] is a real matrix satisfying 00 T = I = T 0. It is an element of the 
rotation group 0(d 2 — 1). Now consider 

*U-< W = t <p'K ® K ® ■ ■ ■ ® \ N ) 

= Trip^l.U^K ®\ 2 ®---® K N (®liU})) 
= TripUiXiM ® U 2 X t2 U\ (8) ... (8) U N X iN U ] N ) 



= J2 Tr (p^ 1 ®x a2 ®---®x aN )olli i oi? a2 ---o 

Oil •••OtN 

= V t o (1) o (2) • o 



/*1*2"'*JV 



where is the k-mode product of a tensor e IR^ 2-1 )*^ 2-1 )*'''*^ 2-1 ) by the 

orthogonal matrix {k) E R^ 2 - 1 )*^ 2 " 1 ) [46,47,49]. Therefore, 

r iN) = xx 0« x 2 o( 2 ) • • • 
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By proposition 3.12 in [46] we get 

I \r {N) 1 1 = \\tw Xl o«x 2 o& . . . x N o w 1 1 = | |rw i | 



Proposition 6 : If a multipartite pure state |^) is subjected to a local measurement 
on the /cth qudit giving outcomes i& with probabilities pi k and leaving residual iV-qudit 
pure state \4>i k ) then the expected entanglement ^2 ik Pi k Eq-{\4>i k )) of residual state is not 
greater then Er(\ip))- 

Eft*Br(W)<W)- (20) 

Proof. Local measurements can be expressed as the tensor product matrix D = 
(g>I)( 2 ) ® • • -<g) D( N ^ on the expanded coherence vector T [50]. The expanded coherence 
vector T is the extended correlation tensor T (defined below) viewed as a vector in the 
real space of appropriate dimension. The extended correlation tensor T is defined by the 
equation 

^ d 2 -i 

p= 2* r hi2-i N K®K®---®^iN ( 21 ) 

where A^ G {/, Ai, A2, • • • , A^-i} are the i^th local SU(d) generators acting on the 
kth qudit (Ao = I) and the real coefficients T ili2 ... iN are the components of the ex- 
tended correlation tensor T. Eq. (6) and Eq.(21) are equivalent with Tqoo-o = 1 > 

/t- _ (1) rr _ ^-{1,2,-M} , _ • • • / n 

AiOO-0 — Su , • • • , ^iii 2 -iM00-0 — 1 i 1 i 2 -i M 5 ' " " ana 1 hi2-i N ~ 1 i l i 2 ..-i N ^ ' ' ' > ? JV 7^ U - 

£)W; fc = 1,2, •••iV are d 2 x d 2 matrices. Without losing generality, we can assume 
the local measurements to be POVMs, in which case £>(*) = diag(l,DW) and the 
(d 2 - 1) x (d 2 - 1) matrix is contractive D^ T D^ < I [50]. The local POVMs 
acting on a iV-qudit state p corresponds to the map p i — > -M(p) given by 

M( P )= y: i^i?®---®^®^®---®^ 

iii2—ijv 

where are the linear, positive, trace preserving operators satisfying Yli k L^L^ = I 
and = 0. The resulting correlation tensor of M.(p) can be written as 

T'W=TWx 1 DUx 2 DW...x N DW 

where is (d 2 - 1) x (d 2 - 1) matrix and D^ T D^ < I. 

Action of POVM on kth qudit corresponds to the map Mk(\ijj) = J2i k Mi k pMj k 
where M lk = I ® ■ ■ ■ ® • • • I , £. L^g' = / and [4 fe) , 4 fe)t ] = with the resulting 
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mixed state £V Vi k \ < t>i k )( < t>i k \, where \<p ik ) is the iV-qudit pure state which results after the 
the outcome i& with probability p ik . 

The average entanglement of this state is 



ik ik 



(AT) 1 1 /^(^-l)\(JV/2) 
)ll ( 2 J 



^ ft J|7^ Xfc D(*)||-(^^)(^) 



(JV/2) 



where, by proposition 3.7 in [46], D^T^(\ ip)) is the fcth matrix unfolding [34] of %^ Xfc 

D^ k \ Therefore, from the definition of the Euclidean norm of a matrix, \\A\ \ = y^Tr(AA^) 
[51] we get 



ik ik 
ik 



= I ITT 



(AO i 



;^-^)W 2 ) = £ r (hA>) 



» ii v 2 

because D^ T D^ < I, and = 1. We have also used the fact that Euclidean norm 

of a tensor equals that of any of its matrix unfoldings. ■ 

Proposition 7 : Let \ip) be an iV-qudit pure state. Let p denote the reduced density 
matrix after tracing out one qudit from the state \ip). Then 

ll^ll^lldl 

with equality only when \ip) = [</>) (g) \x) where \x) is the state of the qudit which is traced 
out. 

Proof. We prove this for a case of qutrit states whose generalization to all qudit 
states is straightforward. Consider a N-qutrit state \ip) = aj6i • • • 6jv) + b\b[ ■ ■ ■ b' N ) + 
c\b'(---b%); |a| 2 + |&| 2 + |c| 2 = l. 
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Here \bi), and \b") are the eigenstates of operating on the ith qutrit. Now 
consider set S of iV-fold tensor products of qutrit operators , namely S = {\ ai ® A Q2 <g> 
• • • <E> A ajv }, «i • • -a N — 1, • • • ,8. 

Choosing , = 8 we get 



(-2Y 

A 8 <g> A 8 <g> • • • <g> A 8 |6i • • • 6jv) = ^]v/r|6i 



C-2) s ' 

A 8 ® A 8 <g> • • • ® A 8 |&i • • • = ^H&i • • • b' N ), 



As®A 8 ®---®A 8 |6' 1 '-..6';) = y^|6 / 1 / ---^), 

where s, s', s" are the number of occurrences of state |3) in |6i • • • 6jv) 5 l^i " " " ^jv) anc ^ 
| b'{ ■ ■ ■ b" N ) respectively. We can replace k of the (N — s), (N — s'), (iV — s")A 8 operators 
acting on the states |1) and |2) occurring in \b± ■ ■ ■ b^), • • • b' N ), \b" ■ ■ ■ b" N ), respectively 
by A 3 to get ( N k s ), { N k s ), ( N ~ k s ) new operators which we call Y k (k — 1, 2, • • • , N — s), 
Yl (k = 1, 2, • • • , TV - s') and F fc " (Jfe = 1, 2, • • • , N - s") respectively. 

We can choose an operator from S, denoted B, such that B\b\ ■ ■ ■ bjq) — i^l&'i • • • b' N ) , 
where 77 is determined by the eigenvalues of A 8 . If B contains q < N occurrences of 
operators Ai, A4, A6, then we can replace k < q of them by A2, A5, A7 operators re- 
spectively. We denote the tensor product operator obtained after k such replacements by 
B k (Bq = B). We have, B k \b\ ■ ■ - b N ) = ±(i) k r]\b' 1 ■ ■ -b' N ). We also choose an operator E 
from S such that E\b\ ■ ■ - b N ) = ±x\b'( ■ ■ -b%), also the operators E k (E = E) (k < q') 
such that E k \b\ ■ ■ ■ b^) = rk{i) k x\b'[ ■ ■ ■ b" N ). Further, we choose an operator D in S such 
that D\b[ ■ ■ ■ b' N ) = ±5\b'[ ■ ■ ■ b%), and also D k (D = D) (k < q") such that D k \b\ ■ ■ ■ b' N ) = 
±(i) k 5\b" ■ ■ ■ b" N ). We can Now calculate various tensor elements, 



(61 • • • &at|A 8 ® • • • ® A 8 |6i ■■■b N ) = 
(b' 1 ---b' N \\ 8 ®---®\ 8 \b' 1 - = 

{b'[---b" N \\z®---®\z\b'[---b" N ) = - 



-2y 



3 N/2 

3^/2 
-2Y" 



3W/2 • 

Replacing (N — s), (N — s'), (N — s")\ 8 by operators by A 3 operators we get 

-2Y 



(61 ■ • -&jv|A 8 <g> ■ ■ ■ <g> A 8 <g> A 3 • • • <g> A 3 1 61 ■••b N ) 
(b[ ■ ■ ■ b' N \X 8 <g> • • • <g> A 8 <g> A 3 • • • <g> A 3 |6; • • • b' N ) 
(b'[ ■ ■ ■ b" N \X 8 <g> ■ ■ ■ ® A 8 ® A 3 • • • ® A3I61 • • • 6^) = 



3^/2 

(-2)" 
3 s '/2 

[-2f 
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The symbols q' and q" have the same definition as that for q. Symbols s' and s" denote 
the number of times the state |3) has occurred in \b[ ■ ■ -b' N ) and \b" ■ ■ - b" N ) respectively. 
The contributions of operators B, B k , E, E k , D, D k are, 

(h • • .ftjvlBl&i • • -b' N ) = ± V = (&i • ■■b' N \B\b 1 ■■■b N ) 
(b' 1 ---b' N \B k \b 1 ---b N ) = ±(i) k r ] 
(b 1 ---b N \B k \b' 1 ---b' N )=±(-i) k r ] 

(&!••• b N \E\b'[ ■ ■ ■ b" N ) = ± X = (K ■ ■ ■ b%\E\h ■■■b N ) 
(b'[---b" N \E k \b l ---b N ) = ±(i) k x 
(b 1 ---b N \E k \b'l---b%) = ±(-i) k X 



also 



(b[ ■ ■ ■ b' N \D\b'[ ■ ■ ■ b" N ) = ±5= (b'[ ■ ■ ■ b" N \B\b\ ■ ■ ■ b' N ) 
(b';---b%\D k \b[---b' N )=±(t) k 5 
(b[---b' N \D k \b';...b%) = ±(-i) k 5 
Finally, we get , for the tensor element, 



3JV 



3 N 

= 7fNM 2 (b l ---b N \\ ai ®---®\ aN \b l ---b N ) + \b\ 2 (b' 1 ---b' N \\ ai ®---®\ aN \b' 1 ---b' N ) 



+\c\ 2 (b'[---b%\\ ai ®---® K N \b'[---bl 



Nl 



+a*b(h ■ ■ ■ b N \ A ai <8> • • • <g> A ajv |&i • • • 6^) + afe* (&i • • • 6^1 A ai <g> • • • <g> A ajv |&i • • • b N ) 
+a*c(6x • • • b N \ A ai (8) ••• (8) A a J&;' • • • 6^) + ac*<&? • • • &^|A ai <g> ■ ■ ■ <g> A Q J&i ■ ■ ■ b N ) 
+Wc{b\ ■ ■ ■ b' N \\ ai (g) ... (g, \ aN \b'[ ■ ■ ■ b" N ) + bc*(b'; ■ ■ • 6';|A Q1 ® • • • <g> A a J&i • • • 6^>] 
The nonzero elements of t ai ... aN are t 88 ... 8 = -^((-2) s |a| 2 + (-2) s '|6| 2 + (-2) s "|c| 2 ), 

^88-83-3 = l a ! 2 + ^ 3 N-k \b\ 2 + { J^l-k l c H> where k, is the number of A 3 s in the element 

A ai ® • • • ® A ajv di = 3, 8. The elements of corresponding to B type of operators are 
t B = ±i]ab* ± i]a*b = ±2i]\a\ \b\cos((fi a — <f) b ), 



t Bk = ±(i) k r)ab* ± (-i) k r)a*b = | 



\k„„& _i_ (_i\k„„*u _ j ±2^|a| \b\cos(<p a - 4> b ) if k is even 

±2?7|a| \b\sin((f) a — (pi,) if k is odd 



We get Ylk=o (2I) elements with cos((fi a — 0&) and X^/Lo Gfc+i) elements with sin(<fi a — 
(fib). If g is odd (for the given stated)) the number of cosines and the number of sines are 
equal. When q is even the number of cosines exceeds by 1. Similarly 
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t E = ±xab* ±xa*b = ±2x\a\ \b\cos{(p a - C ), 

+ - +a\kr* , / _ / ± 2 xM |6|cos(0 a - C ) if fc is even 

t Ek - ±(.) X a6 ± (-.) X a 6 - | ±2x|fl| |6|flin( ^ _ ^ if fc ig odd 

We get Ylt=o (21J e l emen ts with cos(0 a — C ) and YH=o (2^+1) e l emen ts with sin((f) a — 
C ). If g' is odd (for the given state l^)) the number of cosines and the number of sines 
are equal. When q' is even the number of cosines exceeds by 1. 

t D = ±5ab* ± 5a*b = ±2S\a\ \b\cos((f) b - C ), 

t Dk = ±(i) k 5ab* ± {-i)Ha*b = I fj x \ a \ l\ CO f» - *°\ * \ 18 
k w v 1 \ ±25\a\ \b\sm((j)b - C ) if k is odd 

We get Ylt=o ilk) e l emen ts with cos(0 fe — C ) and Ylk=o Gfc+i) e l emen ts with sm(06 — 
C ). If g" is odd (for the given state^)) the number of cosines and the number of sines 
are equal. When q" is even the number of cosines exceeds by 1. 



Finally we get, 

\q-( N )\\2 _ (3^\2 v J_ 
\ 2 \iP) N ~ \2 N ) X 3 N 



((-2)>| 2 + (-2) s >| 2 + (-2) s "|c| 2 )' 



+ e '< * s > e K 1 8 

k=l v 7 fc=l v 7 fe=l v 

+V|a| 2 |6| W(0 a - b ) £ ( 9 ) + V|a| 2 |6| 2 S m 2 (0 a - b ) ]T ( 9 

fc=0 ^ 7 fc=0 ^ 

+4 X 2 |a| 2 |c|W(0 a - C ) £ + 4 X 2 |a| 2 |c| 2 S m 2 (0 a - C ) £ f 9 ' 

fc=0 ^ 7 fc=0 ^ 

+45 2 |6| 2 |c| W(0 6 - C ) ]T (Q + 45 2 |6| 2 |c| 2 W(0 fe - C ) ]T U' x ) 

fc=0 ^ 7 fe=0 ^ 7 



Note that, using |a| 2 + |6| 2 + |c| 2 = 1, it is easy to see that ||7j^|| > 3 N ^ 2 , showing 
that E T > 0. 

Next we consider 

IVXVI = M>i • • • 6^>(6i • • 'M + \b\ 2 \b[ ■ ■ -6^><6i • • -^1 + |c| 2 |6f • • .^(ftf • • -6^| 

+ac* |6x • • • 6jv) (6'/ • • • + a*c|6? • • • b" N ) (h---b N \ 
+bc*\b' 1 ■ ■■b' N ){b'i ■ ■■b'^\ + b*c\b'[ ■ ■■b" N )(b' l ■ ■ -b' N \ 
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and trace out the iVth qudit to get the N — 1 qudit reduced density matrix 
p = \a\ 2 \b x ■ ■ ■b N - 1 )(b l ■ ■ -^-il + \b\ 2 \b[ ■ • •6' JV _ 1 )(6' 1 • • A-il + |c| 2 |&i •••^-iX&i- 

' & jv-iK & Jv!&jv> + a*c|&" • ■■b" N _ 1 ){b 1 ■ ■ -b N -i\{b%\b N ) 

■K-iWnK) + b*c\K ■ ■ ■■b'v-MK) 



■b'^l 



+ab*\b 1 ---b N - 1 ){b' 1 
+ac*\b 1 ■ ••b N - 1 ){b" 

+bc*\b , 1 -.-b' N _ 1 )(b / ; 



Now 



t 

+\b\ 



ai---ajv-i 



= Tr{p\ ai <g> \ a2 


8> • • • A ajv _ 


J = 






-i|A ai O • • • O A ajv _ 1 


61 • • -&jv-l) 


•&JV-l|A Q1 <8>- ■ - ®A 


ajv-i 1^1 ' ' 




1 ) + N 2 ( 


5;'... 


^_i|A Q1 <g)- • - OA ajv 




+a*b(b 1 ■ ■ 


&AT-l|A ai 


® • 


• O A aAr _ 




■■b' N ^)(b N \b' N ) 




+ab*(b[-- 


b^-i A ai 




• O A aiV _ 


Jh- 


■ • b N -i) (b' N \b N ) 




+a*c(b 1 ■ ■ 


&AT-1 |A Ql 




• ® A ajv _ 




• -6jv-i)(6iv|6jv) 




+ac*(b'l-- 


^JV-l|Aai 


® • 






• •&jv-i)(&jv|&jv) 




+b*c(b[ ■ • ■ 


^AT-l|Aai 


® • • 


• <8> A aiV _ 




••^-1><«> 




+bc*(b'l-- 


&JV-1 A Ql 




• <S> A ajv _ 




••6iv-i)(6jv|6jv) 





Wehavefor X— 1 tensor product operators As^Ag®- • -CgiAsl&i • • -6jv-i) = 3 (iv-i)/2 l^i 



1 6l ' • • &JV-.I 



We construct the operators X, X fc , X', X£, X", X^', G, G k , Q, Qk, Z and 
Z k corresponding to Y, Y k , Y', Y k \ Y" , Y k " , B, B k) E, E k) D and D k respectively 
acting on X — 1 qudits. We then get G\b\ ■ ■ - b N ) = ±77 1 6^ • • ■& / A r_ 1 ), G k \b\ ■ ■ ■ 6jv-i) = 
±(i) k r ] '\b' 1 ---b' N _ 1 ) 

Q\h--- b N ) = ± X '\b'l ■ ■ ■ b" N ^), Q k \h ■ ■ ■ ftjv-x) = ±(i) Y|6? • • • 6^_!) 
Z|6i • • = ±8>\b>{ ■ ■ .6^_ 1 >, • • -6^) = ±(i) fe 5K • 



Now, the nonzero elements of X 



(JV-l) 



are U 



[(-2Y\a\ 2 + (-2Y'\b\ 2 + (-2f\c\ 2 }, 



where £, £', £" are defined as s, s', s" for 

t G = ±r]'ab*(b N \b' N ) ±r]'a*b(b' N \b N ) = ±2r t '\a\\b\\(V N \b N )\cos(<f> a -<j> b - a G ), 

t Gk = ±(i) k r ] 'ab*(b N \b' N } ± (-i) k V 'a*b(b' N \b N ) 



±2i]'\a\ \b\\(b' N \bN)\cos((f) a — (fib — oca) if k is even 
±2r)\a\ \b\\{b' N \b N )\sin((fi a - <fi b - a G ) if k is odd 

t Q = ±x'ac*(b N \b%)± X 'a*c(b%\b N ) = ±2 X '\a\ \c\ \ (b%\b N ) \cos{<fi a - <fi c - a Q ), 



t Qk = ±{i) k x!ac*(b N \b" N ) ± (-t) k X 'a*c(b%\b 



N/ 
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±2x'|a| \c\\(b%\b N )\cos((f) a — (j) c — a Q ) if k is even 
±2x|a| \c\\(b%\b N )\sin((f) a - C - a Q ) if k is odd 

t z = ±5'bc*(b' N \b%) ± 5'b*c(b" N \b' N ) = ±25'\b\\c\\(b" N \b' N )\cos{ ( t> b - & - a z ), 
t Zfe = ±(i)^'6c*(^|6^) ± {-i)H'b*c{b" N \b' N ) 



±25'\b\ \c\\(p" N \b' N )\cos(4>b — 4> c — Oiz) if k is even 
±2<J'|6| |c||(6^|6^)|sm(0 6 -0 c -Q; z ) if fc is odd 



Finally we get 

iS°n 2 = (§^) 2 x^[(-2)>i 2 +(-2)nfci 2 +(-2rici 2 ) 2 +Ef- 1 i - s (afe^)Vi 2 rr s ) 



fc=i v 7 fc=i v 



N — 1 — s" 

k 



+4 V ' 2 \a\ 2 \b\ 2 \(b' N \b N )\ 2 cos 2 (<l )a - 0, - a G ) ^ 



fc=0 



m 
2k 



+4r ] ' 2 \a\ 2 \b\ 2 \(b' N \b N )\ 2 sin 2 ( ( j )a - ( f> b 



+ 4X-' 2 |0.| 2 |C| 2 |(( ) ;|V>I' 2 C0S 2 (^ - & - Cg) JJ 



k=0 



m 
2k 



+4 X ' 2 H 2 |c| 2 |(^M| 2 S m 2 (0 a - C - a g ) £ 



fe=0 



2A; + 1 



+45' 2 |6| 2 |c| 2 |(6';|^)| 2 co S 2 (0 fe -<j> c - a z ) £ 



fc=0 



+45' 2 \b\ 2 \c\ 2 \(b%\b' N )\ 2 Si n 2 ( ( p b -<f> c - a z ) ]T 



fc=0 



m 
2k 

to" 
2A; + 1 



where to < g is the number of Ai operators in G, m! < q' is the number of A4 operators in 
Q and to" < q" is the number of A 6 operators in Z. Since 1(^1^) | 2 < 1, |(^jvl^)| 2 < 1 
and |(6jvI^jv)| 2 < 1 we see that 

II^I^H^H 2 
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equality occurring when |6jv) = |&jv) m which case |-0) 
ward, but tedious to elevate is proof for the general case 



10) <8> |6jv) - It is straightfor- 



1,2, 



,8 



ai"-ajv 



Basically we have to keep track of Q V, -B, E 1 , D type of operators, where r is the 
number of terms in the expansion of in order to obtain all nonzero elements of 77 



When iVth particle is traced out, the corresponding elements of Tp N ^ get multiplied by 
the overlap amplitudes, which leads to the required result. The generalization to qudits 
requires action of A operators on the computational basis . Rest of the proof has a straight 
forward generalization. ■ 



(JV) 



Continuity of Er'- We show that for iV-qudit pure states Hd^)^! — |0)(0|)|| — > 



E T {m - E T {\<t>)) 
Proof. 
But 

Therefore 







-10X01)11-0 



2 \4>) I 



> 



\ 2 w 1 



j-(N) _ r (N)\ 



E T (m) - E r {\4>)) 







l T W, 
\ I \4>) I 



I^V) I 



_ iir^i 
- W 1 ^) I 



4.1 Entanglement of multiple copies of a given state 

LU invariance. We show that Ej- for multiple copies of iV-qudit pure state \ip) is LU 
invariant. Consider a system of iV x k qudits in the state \x) = \^) <S> \^) <8> • • • <8> \ip) (k 
copies). It is straightforward to check that [34] 

This implies, in a straightforward way, that 

IM| X > II - lM|v-> II • 

Since by proposition 6 ||7j^|| is LU invariant, so is ll^j^ll- 

Superadditivity : We have to show, for iVqudit states \ip) and |0) that 

E T (m ® |0)) > #r(|V» + ^r(]0))- (23) 

We already know that for \x) — <S> |0) 

IItWii _ iitWii iitWii 
ll^lx) II ~ H J W II \n*) II 
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Thus Eq. (23) gets transformed to 

HT^II IIT (iV) ll - ( d<yd ~ l \ N > IIT (7V) II -i- IIT (iV) ll - n( d ( d ~ l ) \N/2 
W 1 ^) N ll^> III 2 J - ll J |v-> II + N^> II ^ 2 ' 

Putting H^ll = [(fcll)iv/2 + a .j. ^ > this inqua i ity reduces to (^bl))^/2 > 1 - § 
which is satisfied by all d > 2 for all x > 0. 

4.2 Computational considerations 

Computation or experimental determination of Eq- involves (d 2 — 1) N elements of so 
that it increases exponentially with the number of qudits N. However, for many impor- 
tant classes of states, Eq- can be easily computed and increases only polynomially with 
N. We have already computed Eq- for the class of N qutrit GHZ states. For symmetric 
or antisymmetric states is supersymmetric, that is, the value of its elements are 

invariant under any permutation of its indices [34]. This reduces the problem to the com- 
putation of ^U 7 k=1 (N + k) distinct elements of T^ N "> for qutrits, which is a polynomial in 
N[52]. 



5 Extension to mixed states 

The extension of Eq- to mixed states p can be made via the use of the convex roof or (hull) 
construction as was done for the entanglement of formation [16]. We define Er{p) as a 
minimum over all decompositions p = i^o pure states i.e. 

E r (p)= min V ft £ T (|^)). (24) 

The existence and uniqueness of the convex roof for Eq- is guaranteed because it is a 
continuous function on the set of pure states [53] . This entanglement measure is expected 
to satisfy conditions (a), (b) and (c) given in section 4 and is expected to be (d) convex 
under discarding of information, i.e. 

Y^PiMpi) > EriYjPiPi)- (25) 

i i 

The criteria (a)-(d) above are considered to be the minimal set of requirements for 
any entanglement measure so that it is an entanglement monotone [44]. 

Evidently, criteria (a) and (b) are satisfied by Eq-{p) defined via convex roof as it is 
satisfied by Eq- for pure states. Condition (d) follows from the fact that every convex 
hull (roof) is a convex function [54]. We need to prove (c), which is summarized in the 
following proposition. 
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Proposition 8: If a iV-qudit mixed state p is subjected to a local operation on i th 
qudit giving outcomes k with probabilities p k and leaving residual iV qudit mixed state 
Pk, then the expected entanglement 'Y^ k PkEr{Pk) of the residual state is not greater than 
the entanglement Eq-(p) of the original state. 

Y^PkEriPk) < Er(p) 
k 

(If the operation is simply throwing away part of the system, then there will be only one 
value of k, with unit probability.) 

The proof follows from the monotonicity of Er(\ip)) for pure states that is propositions 
5,6 and 7. Bennett et al. prove a version of proposition 8 in [48], which applies to any 
measure satisfying propositions 5,6 and 7. Thus the same proof applies to proposition 8, 
so we skip it. 

Note that any sequence of local operations comprises local operations drawn from the 
set of basic local operations (i)-(iv) stated in section 4, so that proposition 8 applies to 
any such sequence. Thus we can say that expected entanglement of a iV-qudit system, 
measured by 3r(p), does not increase under local operations. 



6 Summary 

In summary we state that, we propose an entanglement measure (Eq. (8)) for iV-qudit 
pure states which passes all the tests expected of a good entanglement measure. More- 
over, this measure is experimentally obtainable, without a prior knowledge of the state of 
the system, by measuring A operators ( generators of SU(d)) which are simply related to 
the angular momentum operators [54]. This measure can be computed efficiently, at lest 
for the state belonging to the symmetric or antisymmetric subspaces. 
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